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Abstract 

It is shown that using the observed extensive air shower (EAS) electron and truncated muon 
size spectra at sea level one can solve the EAS inverse problem - reconstruction of primary 
energy spectra and elemental composition, for not more than 2 kinds of primary nuclei. 
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1 Introduction 



High accuracy of modern EAS experiments in the primary energy region of 10 15 eV [1-4] increased 
a number of publications on the solutions of the EAS inverse problem, which is reconstruction of 
primary nuclei energy spectra based on the detected EAS parameters at observation level [5-9]. 
However, discrepancies in results continue to grow. In most of the cases, it is a result of the 
hidden experimental systematic errors due to uncertainties in the response functions of detectors. 
The uncertainty of A-A interaction model at these energies also contributes to the discrepancies 
in primary nuclei energy spectra results. However, there are certain publications where the EAS 
inverse problem is solved based on erroneous presuppositions, which of course only increase the 
existing discrepancies of data. 

2 EAS inverse problem 

In general, the relation between energy spectra (d?SA/dE) of primary nuclei (A = 1,4,... 56) and 
measured EAS electron and muon size spectra at observation level (AI / AN*^) at a given zenith 
angle 9 is determined by an integral equation 



where the EAS size spectra dW/dN e ^ depend on observation level, A, E, 9 parameters of primary 
nucleus and A — Aait interaction model, dG(N, N*)/dN* is the error function of measurements. 
Equation (1) is a typical ill-posed problem and has an infinite set of solutions for unknown primary 
energy spectra. However, the integral equation (1) turns to a Frcdholm equation if a kind of a 
primary nucleus is defined directly in the experiment (as it successfully does in the balloon and 
satellite measurements [10, 11] where energy spectra for different nuclei are obtained up to 10 15 
cV) or the detected EAS parameters (or combination of parameters) slightly depend on primary 
nuclei [5]. In the last case one evaluate only the all-particle energy spectra ^2 A d?SA/dE from 
equation (1). 

As a result, the meaning of data presented in the [9] and further publications [12-14] where authors 
claim to have got the solution of equation (1) for 4 types of primary nuclei, is not clear. 
Let's prove that using the measured EAS electron and truncated muon size spectra at 3 zenith 
angular intervals, as its done in [9], the Eq. (1) can have a single solution only for primary flux 
which consists of not more than 2 kinds of nuclei. 

Let the EAS size spectra at observation level KASCADE (t = 1020 g/cm 2 [9]) be described by 
log-Gaussian distributions with mean < N e ^(E, A, 9) > and variance er 2 ^(E , A, 9) . It is known 
that this assumption is well performed at atmosphere depth t > 700 g/cm 2 , primary energies 
E > 10 5 GeV and zenith angles 8 < 35°. 




(1) 



The kernel function of equation (1) is determined as 




(2) 



Let also the measurements and further evaluations of EAS electron and truncated muon sizes be 
carried out without errors (dG/dN* = S(N — N*)) and integral Eq. (1) include only statistical 
uncertainties. 

Then a set of Eq. (1) transforms into the following: 
AI(6) ^ fdZ A dW(E,A,0) 

Let's determine the parameters of distribution functions dW/dN e ^ by the following known em- 
pirical expressions: 

<N>^a( 1 £ v ) i A' c „ d e (4) 

where the values of corresponding approximation parameters (a, ... d, a, ... p) are presented in 
Tables 1,2 and obtained by CORSIKA6016(NKG) EAS simulation code [15] at QGSJET interaction 
model [16]. 



a 


b 


c 


d 


e 0.0041±0.0004 


1.22±0.005 


-0.22±0.01 


6.8±0.2 


H 0.0079±0.0003 


0.933±0.001 


0.07±0.002 


1.87±0.05 



Table 1: Approximation parameters of average EAS electron size (N e ) and EAS truncated muon 
size (iV M ) in empirical formula (4). 



a 


S 


e 


P 


e 0.79±0.01 


0.25±0.003 


1.82±0.05 


2.1±0.1 


0.283±0.003 


0.277±0.003 


1.12±0.05 


0.26±0.1 



Table 2: Approximation parameters of a e ^ in empirical formula (4). 



The accuracy of approximation (4) for average EAS electron size < N e (E,A,6) > is less than 
10% at 3 • 10 5 < E < 3 • 10 8 GeV, A = 1, ... 56, 6 < 32° and observation level 1020 g/cm 2 . The 
corresponding accuracies of cr e , M and average EAS truncated muon size AT M are less than 1-2%. 
Changing the variables of kernel functions of Eq. (3) from N e and to x e and respectively 
according to 

we obtain the following set of Fredholm integral equations 



AJ(fl) 
AA e ^ 



- j f eitl (E,A)F e ^(E,0)dE (7) 



where 

f eA E,A)=J2^i AS ^ C ^ ( 8 ) 

A 

and kernel functions 

independent of kind of primary nuclei. 

Examples of distribution functions dW/dx for EAS electron and truncated muon size spectra are 
presented Fig. 1,2 respectively. These data were obtained by CORSIKA code for primary proton 




Figure 1: EAS electron size spectra for primary proton (empty symbols) and iron (filled symbols) 
nuclei at energy E = 3.2 PeV, zenith angle 9 < 18° and observation level 1020 g/cm 2 . The lines 
are the corresponding log-Gaussian distributions. 

(A = 1, empty symbols) and iron (A = 56, filled symbols) nuclei at energy E = 3.2 • 10 6 GeV and 
zenith angle < 18°. The lines in Fig. 1,2 correspond to log-Gaussian distributions which were 
counted based on known values of mean < x >= 1 and variances a^ ^(E, 9). 

It is seen, that the distribution functions slightly depend on a kind of a primary nucleus, especially 
the right-hand sides of distributions. It is an important fact, because the contribution of the left- 
hand sides of distributions is negligible small at a priori steep primary energy spectra (d^SA/dE ~ 
E~ 3 ). 

Let the functions fe*' k \Ei) and fjf' k \Ei) be the solutions of the set of Fredholm equations (7) at 
i = 1, . . . m values of primary energies Ei, j = 1, . . . n kinds of primary nuclei and k=l,2,3 zenith 
angles. Then the values of unknown energy spectra (d^AjdEi) at a given set of Ei and different 
primary nuclei (Aj = Ai, . . . A n ) are determined by the set of linear equations 

<&\E i )>e= E^-^L,, : i = l,-m (10) 

A=Ai * 



^4. — A n 

<f { j ) {E i )> e = ]T AS 

A=A X 

rU)t 



dE 



E=E, 



i = 1, . . . m 



(11) 



where < > e = § £ fc /i?£ fc) (£i) 

Evidently, the single solutions of linear equations (9,10) occur only at ran < 2m, where ran is a 
number of unknown (d^SAj/dEi) and Ira is a number of equations (10,11). 
Therefore the condition of existence of solutions is n < 2 and the energy spectra [9] obtained 




Figure 2: The same as Fig.l for EAS truncated muon size. 

based on KASCADE EAS electron and truncated muon size spectra for 4 (A = 1, 4, 16, 56) primary 
nuclei no physical meaning. At the same time the all-particle spectrum obtained in [9] have to be 
approximately right due to the power index in Eq. (11) 5^ — <C 1. 

It should be noted that unreliability of solutions [9] is also shown in [17] based on correlation 
analysis between the observable electron and truncated muon size spectra and unknown energy 
spectra of primary nuclei. 



3 Conclusion 

Based on equation (1), the investigation of the EAS inverse problem, which is the reconstruction 
of energy spectra of primary nuclei by the observable EAS electron and muon size spectra at 
observation level makes sense only using a priori given functions for unknown primary energy 
spectra with given unknown spectral parameters (so called parameterization of equation (1)) as it 
is done in [8, 18, 19, 20]. 
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